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A gravitational wave traversing the line of sight to a distant source produces a frequency shift 

s : 

which contributes to redshift space distortion. As a consequence, gravitational waves are imprinted 
as density fluctuations in redshift space. The gravitational wave contribution to the redshift space 
power spectrum has a different u dependence as compared to the dominant contribution from 

U: 

peculiar velocities. This, in principle, allows the two signals to be separated. The prospect of 
a detection is most favourable at the highest observable redshift z. Observations of redshifted 

o ■ 

21-cm radiation from neutral hydrogen (HI) hold the possibility of probing very high redshifts. 
We consider the possibility of detecting primordial gravitational waves using the redshift space HI 

, power spectrum. However, we find that the gravitational wave signal, though present, will not be 

>' 

detectable on super-horizon scales because of cosmic variance and on sub-horizon scales where the 
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^sO ■ signal is highly suppressed. 
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I. INTRODUCTION 



Primordial gravitational waves are a robust prediction of inflation [1, 2]. These stochastic 
tensor perturbations are generated by the same mechanism as the matter density fluctua- 
tions, the ratio of the tensor perturbations to scalar perturbations being quantified through 
the tensor-to-scalar ratio r. Detecting the stochastic gravitational wave background is of 
considerable interest in cosmology since it carries valuable information about the very early 
universe. The cosmological background of gravitational wave has its signature imprinted on 
the CMBR temperature [3] and polarization [4] anisotropy maps. Current CMBR observa- 
tions (WMAP-5 Year data) impose an upper bound (r < 0.43) which is further tightened 
(r < 0.22) if combined CMBR, BAO and SN data is used [5]. Detecting the gravitational 
wave background is one of the important aims of upcoming PLANCK [6] mission and future 
polarization based experiments like CMBPOL [7]. 

A gravitational wave traversing the line of sight to a distant source will contribute to its 
redshift in addition to that caused by Hubble expansion and its peculiar velocity. This will 
produce a redshift space distortion in a manner similar to that caused by peculiar velocities 
[8]. The effect arises due to the fact that distances are inferred from the spectroscopically 
measured redshifts. As a consequence, a gravitational wave will manifest itself as a density 
fluctuation in redshift space. In this paper we propose this as a possible technique to detect 
the primordial gravitational wave background. 

While one could consider the possibility of detecting this at low redshifts (z ~ 1) us- 
ing galaxy and quasar redshift surveys, we shall show that the prospects are much more 
favourable if the redshift is pushed to a value as high as possible. 

Observations of redshifted 21-cm radiation from neutral hydrogen (HI) can be used to 
measure the power spectrum of density fluctuations at very high redshifts extending all the 
way to the Dark Ages (30 < z < 200) [9]. Redshift space distortions make an important 
contribution to this signal [10]. We investigate the possibility of using this to detect pri- 
mordial gravitational waves. We note that the imprint of gravitational waves on the 21-cm 
signal from the Dark Ages has also been considered in an ealier work [11]. 
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II. FORMULATION 



The radial component of peculiar velocity introduces a redshift z v = v/c in excess of the 
cosmological redshift which arises due to the expansion of the universe. This distorts our 
view of the matter distribution in the three dimensional redshift space, where the radial 
distance is inferred from the measured redshift. As a consequence the density contrast 
5 = Sp/p measured in redshift space 5 s is different from the actual density contrast 5 r , and 
[12] 

where H(a) is the Hubble parameter and x the co moving distance to the source. We see that 
any coherent velocity pattern (in-fall or outflow) manifests itself as a density fluctuation in 
redshift space. This takes a particularly convenient form in Fourier space if we assume that 
the peculiar velocities are produced by the density fluctuations 5 r . We then have 

A*(k) = (l + V)A^(k) (2) 

where A s and A r are the Fourier transforms of 5 s and S r respectively, / = dlnD/dlna m 
f2^ 6 , D being the growing mode of density perturbations and p — h • k/fc is the cosine of the 
angle between the line of sight n and the wave vector k. It follows that the power spectrum 
of density fluctuations in redshift space P s {k) 1 is related to its real space counterpart P r (k) 
as [8] 

P s (k) = (1 + fp 2 ) 2 P r (k) (3) 
A gravitational wave h a b(x,rj) which is a tensor metric perturbation 

ds 2 = a 2 [c 2 drf - {5 ab + h ab )dx a dx b ] (4) 

makes an additional contribution [13] 



1 n a 

z h = ~n a n b / h' ab (x,rj)dr] (5) 

to the redshift along the line of sight of the unit vector n. Here prime denotes a partial 
derivative with respect to 77, r\ e and 770 respectively refer to the photon being emitted and 
the present epoch when the photon is observed, and x = n(i]o — rj) is the photon's spatial 
trajectory. Considering Zh the gravitational wave contribution to the redshift, we have an 
additional contribution 

6h ~ aHdx {b) 



to 5 s the density contrast in redshift space (eq. 1). Simplifying this using x = c(r)o — i] e ) we 
have 

« = ^ (7) 
We consider the primordial gravitational waves which we expand in Fourier modes as 

/d 3 k 
h ab (k,rj)e lk - x j^ (8) 

and decompose h a b(k,r]) in terms of the two polarization tensors e^ b and e* b as [14] 

MM) = h(k, V ) [e>+(k) + e>><(k)] V^gggl . ( 9 ) 

Here h(k,r)) quantifies the temporal evolution, and h(k,r)) = 3ji(kcr]) / (kcrj) in a matter 
dominated universe, ji being the spherical Bessel function of order unity. The polarization 
tensors are normalized to e^ b e +ab = e* b e xab = 2 and e^ b e xab = 0, Ph(k) is the primordial 
gravitational wave power spectrum [5] and a x (k), a + (k) are Gaussian random variables such 
that 

(K b (k, V )h ab (k', V )) = (2nf5%k-k')h\k, V )P h (k) (10) 

Let us first consider a single Fourier mode of gravitational wave with k along the z 
direction, and represent the line of sight as 

n = sin #(cos (f>i + sin j) + cos 9 k . (11) 

We can than express eq. (7) as 

7 ' 

A s (k,??) = — -sin 2 ^ [cos20a + (k) + sin20a x (k)l J(2i:fP h {k) . (12) 
AaH 

This can be equivalently interpreted with n fixed and the direction of k varying. We use 
this to calculate -P^(k) the gravitational wave contribution to the power spectrum of density 
fluctuations in redshift space 

n 2 



P, s (k)=sin 4 # 



h 



P h {k) (13) 



AaH 

Thus the total power spectrum of density fluctuations in redshift space is 

P\k) = (1 + f^fP r {k) + (1 - ^ 2 fP r h (k) (14) 
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FIG. 1: This shows the ratio r/r at different z. This is predicted to have a constant value ~ 0.16 
on super-horizon scales in the O m = 0.3 LCDM model considered here. 



where P[{k) refers to the terms in { } in eq. (13). Here P r (k) and P[{k) are respectively the 
matter and gravitational wave contributions to the power spectrum of density fluctuations 
in redshift space. Both P r (k) and P[{k) are to be evaluated at the epoch corresponding to 
the redshift under observation. 

The contributions from P r (k) and P[{k) have different fi dependence. This, in principle, 
can be used to separately estimate the gravitational wave and the matter contributions from 
the observed redshift space power spectrum. While the matter contribution is maximum 
when k and n are parallel, the gravitational wave contribution peaks when the two are 
mutually perpendicular. 



III. RESULTS 



We use f = P[(k)/P r (k) to quantify the ratio of tensor perturbations to scalar pertur- 
bations in the redshift space power spectrum. Assuming n s = 1, C 1, the value of f is 
constant on super-horizon scales (kcrj C 1). This value is f — r/4 if fl m = 1, and somewhat 
smaller (Figure 1) with f = 0.16r for fl m = 0.3 in the LCDM model. Gravitational waves 
decay inside the horizon whereas matter perturbations grow on these scales. The ratio f(k) 
is oscillatory and is severely suppressed on sub-horizon scales [kcrj <C 1). 
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The prospect of detecting the gravitational wave signal is most favourable on super- 
horizon scales (k < kn = (cr/)' 1 ). The k range amenable for such observations (Figure 1) 
increases with redshift z (smaller horizon ci]). Observations of redshifted 21-cm radiation 
hold the potential of measuring the redshift space power spectrum in the z range (30 — 
200) [9, 10], where the pre-reionization HI signal will be seen in absorption against the 
CMBR. Gravitational waves will make a ~ r x 16% contribution to the HI signal on scales 
k < /ch- 

IV. FEASIBILITY OF DETECTION 

The cosmological HI signal will be buried in foregrounds [15-19] which are expected to be 
orders of magnitude larger than the signal. The foregrounds are continuum sources whose 
spectra are expected to be correlated over large frequency separations, whereas the HI signal, 
a line emission, is expected to be uncorrelated beyond a frequency separation. While this, 
in principle, can be used to separate the HI signal from the foregrounds, it should be noted 
that the frequency separation beyond which the HI signal becomes uncorrelated increases 
with z and angular scale. This is a potential problem for the detection of the gravitational 
wave signal. In the subsequent discussion we have assumed that the foregrounds have been 
removed from the HI signal. 

The distinctly different /i dependence of the scalar and gravitational wave components 
of the redshift space power spectrum can in principle be used to separate the two signals. 
Expressing the fx dependence [20] as P s (k ) /i) = Po(k) + P2(k)fi 2 + P 4 (A;)/i 4 , the gravitational 
wave component can be estimated using P[(k) = [Po(k) — P 2 (k)]/2. For a cosmic variance 
limited experiment, the error in P 2 (k) and P Q (k) would be SP(k)/P(k) ~ l/y/N^jk) [17, 21- 
24], where N(k) denotes the number of k modes within the comoving volume of the survey. 
Thus N(k) > r~ 2 ~ 10 4 modes would be needed for a detection of the gravitational wave 
signal. 

The number of modes with a comoving wave number between k and k + dk is dN{k) = 
k 2 dkV/(2ir) 2 , where V is the comoving survey volume. Assuming a survey between z = 20 
to z — 200, and using a k bin dk = k/10, we have N(k) = 10 for k = ku ~ 0.002 Mpc -1 . 

It is, in principle, possible to carry out HI observations in the entire z range z = to 
z = 200 [15] and thereby increase the volume. Of the entire survey volume Vo, for a mode k 
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only a volume V(k) = Vo — (47r/3)(c7]o — A:^ 1 ) 3 where the mode is super-horizon contributes 
to the signal. Further, the largest mode k max is the one that entered the horizon at z — 200, 
and the smallest mode k min has wavelength comparable to the radius of the survey volume. 
We then have, assuming a full sky survey, 



which gives N ~ 100. The number of independent modes is too small for a measurement 
at a level of precision that will allow the gravitational wave component to be detected. In 
conclusion, we note that the gravitational wave signal, though present, will not be detectable 
on super-horizon scales because of cosmic variance and on sub-horizon scales where the signal 
is highly suppressed. 
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